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Abstract 

We present a renormalization group theory for the onset of Ising-nematic order in a Fermi liquid in 
two spatial dimensions. This is a quantum phase transition, driven by electron interactions, which 
spontaneously reduces the point-group symmetry from square to rectangular. The critical point is 
described by an infinite set of 2+1 dimensional local field theories, labeled by points on the Fermi 
surface. Each field theory contains a real scalar field representing the Ising order parameter, and 
fermionic fields representing a time-reversed pair of patches on the Fermi surface. We demonstrate 
that the field theories obey compatibility constraints required by our redundant representation of 
the underlying degrees of freedom. Scaling forms for the response functions are proposed, and 
supported by computations up to three loops. Extensions of our results to other transitions of 
two-dimensional Fermi liquids with broken point-group and/or time-reversal symmetry are noted. 
Our results extend also to the problem of a Fermi surface coupled to a U(l) gauge field. 
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I. INTRODUCTION 



A number of recent experiments 1-4 have noted the presence of Ising-nematic order in 
the enigmatic normal state of the cuprate superconductors. This order is associated with 
electronic correlations which spontaneously break the square lattice symmetry to that of a 
rectangular lattice: i.e. the symmetry of 90° rotations is lost, and the x and y directions 
become inequivalent. This broken symmetry is associated with an Ising order parameter, 
which we will represent below by a real scalar field (p. 

Of particular interest are recent experiments on the anisotropy of the Nernst signal 4 in 
YBa 2 Cu30 3/ , which indicate that the Ising-nematic order has its onset at the temperature 
T = T*, which also marks the boundary between the 'pseudogap' region and the 'strange 
metal'. These results call for the theory of the quantum phase transition involving Ising- 
nematic ordering in a Fermi liquid metal. Such a quantum critical point would play an 
important role in the theory of the strange metal. The metallic Ising-nematic critical point 
is also of importance in experiments 5 on S^R^Oy, where the observations of resistance 
anisotropies have demonstrated spontaneous Ising-nematic ordering. Finally, there are clear 
indications of Ising-nematic order driven by electron correlations in the pnictides. 6-9 

One approach to the Ising-nematic ordering is to take a liquid-crystalline perspective 10 , 
and view it among a class of phases with broken square lattice symmetry 11-14 . Ising nematic 
phases are also a generic feature of frustrated and doped antiferromagnets, because the 
Ising-nematic order survives after antiferromagnetism (at wavevectors ^ {^,^)) has been 
disrupted by thermal 15 ' 16 or quantum 17 ' 18 fluctuations. 

A complementary point of view 19-31 is to start from the Fermi liquid with perfect square 
lattice symmetry and look for the Pomeranchuk instability of Landau's Fermi liquid theory 
in the angular momentum £ = 2 channel. Almost all of these works rely on the perspective 
of Hertz 32 , in which the electrons are integrated out to yield a Landau-damped effective 
action for the scalar order parameter 0; the low energy particle-hole excitations near the 
Fermi surface lead to long-range interactions in the action for 0. However, this procedure of 
successive integration of fermionic and then bosonic degrees of freedom is clearly dangerous. 
A systematic renormalization group analysis requires that all excitations at a given energy 
scale be treated together. Consequently, a complete scaling analysis of the Ising nematic 
critical point is lacking: such an analysis should be based on a local field theory, and provide 
a scheme for computing the scaling dimensions of all perturbations of the critical point. 

We can also consider the onset of Ising-nematic order in a superconductor, rather than 
in a Fermi liquid. In a s-wave superconductor, the fermionic excitations are fully gapped, 
and so the theory for <fi has no long-range interactions: consequently the transition is in 
the universality class of the 2+1 dimensional pure Ising model. A ci-wave superconductor 
does have gapless fermionic excitations at special 'nodal points' in the Brillouin zone, and 
these nodal fermions do modify the universality of the transition away from pure Ising 33 ' 34 . 
A fairly complete understanding of the Ising-nematic transition in ci-wave superconductors 
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has been reached in recent work 35 ' 36 using a large- N expansion, where N is the number of 
fermion components. 

This paper provides a scaling theory of the Ising-nematic quantum critical point in two- 
dimensional metals, satisfying the requirements stated above. Our theory builds upon the 
work in the <i-wave superconductor 35 ' 36 , and also on advances by Polchinski 37 , Altshuler, 
Ioffe, and Millis 38 , and Sung-Sik Lee 39,40 on a closely-related problem: the dynamics of a 
Fermi surface with the fermions coupled minimally to a U(l) gauge field. 

We focus on a pair of time-reversed patches on the Fermi surface and describe their 
vicinity by a local 2+1 dimensional field theory. In principle, there are separate critical 
theories for each pair of time-reversed points on the Fermi surface, as is also the case in 
the Fermi surface 'bosonization' methods. 27 ' 28 ' 41 ~ 45 However, a key difference from the latter 
methods is that each Fermi surface point is associated with a 2+1 dimensional theory, and 
not a 1+1 a dimensional theory. This means that there is a redundancy in our description, 
and sowing the theories together is not trivial: we show in Section IV A how this is done in 
a consistent manner. 

Apart from their application to the Ising-nematic transition of interest, simple extensions 
of our results apply also to the U(l) gauge field case, and to other symmetry breaking 
transitions in Fermi liquids involving order parameters which carry momentum Q = 0. We 
will describe these cases in Section II below, and briefly indicate the needed extensions in 
the body of the paper. 

Transitions with order parameters which carry momentum Q ^ lead to different field 
theories, which will be described in a subsequent paper. 46 

After a discussion of the one loop results in Section III, we present our main scaling 
analysis in Section IV. This includes a discussion of Ward identities which strongly constrain 
the structure of renormalization group flow. Finally, explicit three loop computations appear 
in Section V and Appendix B. 

II. THE MODEL 

We consider quantum phase transitions in metals of electrons c a [a =t>4-)> involving an 
onset of a real order parameter 4>{x) at wave-vector Q — 0. The order parameter is taken 
to have the same transformation properties under lattice symmetries and time reversal as, 

9 ka 

For definiteness, we consider a system on a square lattice. Then, <p can describe the following 
patterns of symmetry breaking: 

I. Breaking of the point-group symmetry with = d^, and d^ a = d_^ a . In these cases 
dr has either d x 2_ y 2, d xy , or g-wave symmetry. The Ising-nematic transition of most 
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interest to us here corresponds to the d x 2_ y 2 or d xy cases. These cases all belong to 
one-dimensional representations of the square lattice point group, and we will argue 
that these transitions are all in the same universality class. 

2. Breaking of time-reversal and point-group symmetry with = d^ and d^ c = —d_^ a . 
In this case d^ transforms under the two-dimensional p-wave representation, and so 
requires a two component order parameter = (<p x ,<j) y ). We will not consider the 
two-component case explicitly, but our results have an immediate generalization to 
this transition. This case corresponds to the "circulating current" order parameters 
proposed by Simon and Varma 47 , as was argued in Refs. 33 and 48. 

3. Breaking of spin-inversion symmetry with d^ = —d^,. In this case, dj: can have 
either s-wave symmetry (Ising ferromagnet), <i-wave symmetry (Ising spin-nematic) 
or g-wave symmetry. Unlike transitions i) and ii), which respect the full SU(2) spin 
rotation symmetry, in the present case we assume this symmetry is explicitly broken 
to a U(l) "easy axis" subgroup. 

Notice that in all cases, there is a Z 2 symmetry (either tt/2 rotation, reflection or time- 
reversal) under which <f> — y — <f>. 

Apart from the above symmetry breaking cases, we will also consider the problem of a 
Fermi surface minimally coupled to a U(l) gauge field 37-40 ' 49-58 . This case is similar to case 
2 above, as we describe below Eq. (2.4). Such models arise in theories 57 ' 58 of certain U(l) 
spin liquid phases in which c a describe the fermionic spinons. We will therefore refer to this 
model as the "spin-liquid" case below. The same theory also describes 59-61 "algebraic charge 
liquids" in which case the c a are spinless, charge — e fermions, and a represents the charge 
of the fermion under the emergent U(l) gauge field; we will not refer to this case explicitly 
below. 

Given the order parameter in Eq. (2.1), we may write down an effective spacetime La- 
grangian describing the interactions of the order parameter <fi with the fermions as, 

L = ct(d T + e(-iV))c a - 0{x)<f>(x) + ^(V0) 2 + ^0 2 (2.2) 

Here, we have added by hand a gradient term and a mass for the bosonic mode <fi. Such 
terms will be generated automatically after integrating out the high-energy fermions. The 
absence of higher order terms in <fi and gradients of will be justified below. 

The Lagrangian L in Eq. (2.2) is not yet in a form suitable for our analysis of quantum 
criticality. The main point is that the fermion spectrum e(k) has zeros along the entire 
Fermi surface of large momenta k: so, as is well known, we are not in a position to make 
a low momentum expansion needed for a field theory. One strategy is to use the Hertz 
approach 32 of integrating out all the c fermions to obtain a non-local effective action for the 
order parameter <fi. The latter is singular only at small momenta q and u, and so it is then 
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at least permissible to make a low momentum and frequency expansion. However, the terms 
in the effective for turn out to be highly singular as q — > (see Ref. 26 and Appendix A). 
Moreover, in d = 2, the strength of the singularity increases with increasing powers of <p in 
the effective action. The situation now seems hopeless, but progress becomes possible after a 
key observation: the leading singularities in the <p effective action appear only when all the 
fields have their momenta nearly collinear to each other, as is explained in Appendix A, and 
as will become clear from the structure of our analysis below (by nearly collinear we mean 
that the angle 9 between the momenta is of order 9 ~ \q\/kp). In other words, if we are 
interested only in leading critical behavior, fields with non-collinear momenta effectively 
decouple from each other. The couplings between <p fields with non-collinear momenta are 
then irrelevant corrections to the critical theory. The argument supporting this statement 
is presented in Appendix A. More generally, consider an n-point function 

(0(g 1 )0(g 2 )0(g 3 ) • ■■4>(q n )) ■ 

In a Gaussian theory for 0, which is the claim of Hertz 32 , such a correlator would decouple 
into products over pairs of momenta which sum to zero. However, such a decoupling is too 
drastic: rather, the decoupling is only over sets of momenta which are collinear with each 
other, so that the leading critical singularity of the above correlator takes the form 

II WalMfa) ■■■}■ 
a 

Here all the momenta q a i in a group Q a are collinear to each other, while being non-collinear 
to momenta in groups Qb with b ^ a. We can therefore limit ourselves to <fi fields with 
momenta along a fixed direction q. We will now argue that for each such direction q, there 
is a sensible and powerful continuum limit of Eq. (2.2). 



It is now clear that we may restrict our search for a field theory to that describing the 
singularities in the <p correlations for a single group of collinear momenta Q a . So let us pick 
a direction q for <fi. It is believed that a bosonic mode with momentum q interacts most 
strongly with the patches of the Fermi-surface to which it is tangent 37-40 . Assuming that 
only a single Fermi surface is present, for each q there will be two such points with opposite 
Fermi-momenta ko and —ho, see Fig. 1. We will denote fermions at these momenta as ip + 
and 

4>+v(k) = , V-cr(fc) = C-k 0+ k, a - (2-3) 

We choose coordinate vectors x and y to be respectively perpendicular and parallel to q. 
Then, expanding the fermion energy near k Q and — k , the needed, low energy, continuum 
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FIG. 1. The shaded region represents the occupied states inside a Fermi surface. Fluctuations of 
the order parameter <f> at wavevectors parallel to q couple most strongly to fermions near the Fermi 
surface points ±feo- These fermions are denoted ip±. 

Lagrangian becomes 



Here vf and m are the Fermi velocity and the band mass at ko, while d± tT = d±k oa , and we 
have added a subscript k to L emphasize that this is the Lagrangian for the patch near 



We should emphasize here that all the fields in Eq. (2.4) are 2+1 dimensional quantum 
fields, with full dependence upon x, y, and r i.e. the fields are (p(x, y, r) and ip± a (x, y, r). In 
principle, we should also add a term (d x <f)) 2 to Eq. (2.4); however, we omit it at the outset 
because it will later be seen to be irrelevant near criticality. Further, because of this full 
dependence on x, and y, the fermion fields ip± a describe an extended patch of the Fermi 
surface near the points ±k , and not just the two points ±k . We place some finite cutoff 
A on the size of this patch, and will be interested in the scaling behavior at momenta much 
smaller than this cutoff. 

We now discuss the structure of the couplings d± a in Eq. (2.4). For the transitions in s, d 
and g channels in case 1 above d +a = d_ a by inversion symmetry, and d± a is a independent. 
For case 2, we have d +a = —d- a and also a independent, although the fermions now couple 
to a projection of the two component order parameter • d, while the bosonic gradient term 
generally involves both components of the order parameter. The spin liquid case also has 
d +a = —d_ a and a independent, and cj) is associated with the transverse component of the 




(2.4) 



±k . 
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spatial gauge field in the Coulomb gauge 37 40 ; moreover the spin-liquid has r = by gauge 
invariance. Finally, the Ising ferromagnet case 3 has d +a = gL ct and d±^ = —d±±. 

We note that for transitions in non-zero angular momentum channels, the coupling d 
vanishes along certain axes in the Brillouin zone. The intersections of these axes with the 
Fermi surface are known as cold-spots, as the fermion coupling to the order parameter at 
these points involves additional derivatives and is much weaker. The scaling theory that 
follows only describes the Fermi surface away from cold spots. 

It is convenient to rescale coordinates and fields in (2.4), x = (2mvp)~ l x, ijj = v^ 1 ^ 2 ^, 
(j) = 2m \d\ < f ) - We drop the tildes in what follows. Then, 

L = i\)\ a (r)d T - id x - dfjij; +(T + [r\d T + id x - d 2 ^- a 

- A +CT M> +CT - A_ CT 0^-a + 7^(dy<P) 2 + \tf (2-5) 

with e 2 = 2md 2 /vp, r = ro/(2md 2 ), rj = 2m, and X sa = d scr /\d\, and we will henceforth drop 
the subscript ko on L. We note that as usual, the relation between the parameters of the 
effective theory and the original model should not be taken literally. Rather, in the critical 
regime, we have ro — r$ c = Z r (r — r c ), where r c and tq c denote the critical points of the 
effective theory and the microscopic theory respectively. Moreover, the original fields and 
the fields defined in each patch of the Fermi surface are related by, 

4>(q, u) ~ Z^ /2 K<j) pa t ch {Kq x , q y , u), ip(q, u) ~ Z^ 2 Ktp patch {Kq x , q y , u) (2.6) 

Note that the "metric factors" K, Z r , e 2 , Z^ are generally dependent on the direction 
of the boson momentum q and the cut-off of the low-energy theory A. 

For brevity, we will only present explicit calculations for the case that does not involve 
spin (Ising- nematic transition and spin-liquid); the extension of the results to the Ising fer- 
romagnet case will be noted. Moreover, we extend the number of spin components (flavours) 
to iV from the physical value N = 2 with the view towards performing a large- iV expansion. 
For this purpose, it is convenient to rescale e 2 and r, yielding our Lagrangian in its final 
form 




Here and below we suppress the flavour index. To reiterate, the Ising-nematic case has 
A + = A_ and the spin-liquid case (i.e. Fermi surface coupled to U(l) gauge field) has 
A, = -A_. 
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III. ONE LOOP PROPAGATORS 



To gain some insight into the low energy properties of the theory (2.7), it is useful to 
compute the one loop boson and fermion self-energies. 




FIG. 2. One loop contributions to the (a) boson, and (a) fermion self-energies. 



The one-loop boson polarization in Fig. 2 a) is given by, 

M<l)=N J ^G°(l)G° s (l + q) (3.1) 
We first evaluate this diagram with a bare fermion propagator, 

-i7]k T + sk x + fc* 

The resulting polarization function takes on a characteristic Landau-damped form, 

J {2ir) 2 -vqq T + 2q y l y + q x + q l y 

= ^f^L til , + {q^-q) = c b N^ c b = ^. (3.3) 

2ir J 2ix -iT]q T + 2q y l y + q x + q l y \q y \ An 

Note that r\ has dropped out of the final result. We are interested above only in the singular 
contribution to Ho, and this is insensitive to orders of integration: so unlike the conventional 
order, we have integrated over l x before l T . We include the RPA polarization bubble (3.3) 
into the bosonic propagator to obtain 



- + %L 

NV b \ qy \ + e> 
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Note that the q 2 term is not renormalized by the polarization contribution at this order, and 
the bare co-efficient represents the phenomenological contribution of higher energy modes. 

The one- loop correction to the fermion propagator is given by Fig. 2 b). For simplicity, we 
work at the critical point and set r = 0. Then, the fermion self-energy assumes a non-Fermi 
liquid form 

E.(*) = -/ ^D(l)G° s (k-l) 

__J_ f dlrdly ( jy , 5.V 1 x (k _ r ) 

= -fsgn(MIM 2/3 , =/ = ^(y • (^) 

Note, again, that 77 has dropped out of the result. Incorporating this correction into the 
fermion propagator, 

G s (k) = (- l jL sgn (k T )\k T \ 2 / 3 + sk x + k*\ (3.6) 

Here we have dropped the bare fermion time derivative term proportional to rj, which is 
irrelevant at low energies compared to the dynamically induced self-energy (3.5). 

As is well known, 37 the one-loop expressions (3.3), (3.5) actually satisfy the Eliashberg- 
like equations, in which the lines of Fig. 2 become self-consistent propagators. In what 
follows, we will use these self-consistent propagators (3.4), (3.6) in our calculations and drop 
self-energy corrections like those in Fig. 2. 

IV. SCALING AND RENORMALIZATION 

As has been argued by a numer of authors 37-40 , a useful starting point for the renormal- 
ization group analysis of the theory (2.7) is obtained by using the scaling, 

kg y s k% j y sky j uj y s uj ; 
ip(x, y, r) -»■ s 2 ip(s 2 x, sy, s 3 r), <j>(x, y, r) -> s 2 (f)(s 2 x, sy, s 3 r) (4.1) 

This scaling is suggested by the one-loop calculation of fermion and boson propagators in 
Eqs. (3.4), (3.6). The bare fermion time derivative term ijj^d T ijj is irrelevant under this 
scaling, and so we will take the limit rj — > + . Note that neither of the one loop corrections 
Eqs. (3.3), (3.5) depend upon 77. 

Alternatively, note that the scaling of time in (4.1) could also have been derived by de- 
manding that the 'Yukawa coupling' X s be invariant. This avoids the somewhat unnatural 
appeal to the one-loop self-energy to set bare scaling dimensions, and yields all the scaling 
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dimensions in (4.1) by a simple rescaling of the bare Lagrangian L in Eq. (2.7). Of course, 
once we have set \ s to be invariant, then the coupling 77 becomes irrelevant. These fea- 
tures of the scaling analysis are shared by the theory of the nematic transition in <i-wave 
superconductors in Ref. 36. 

Note also the different scaling of spatial momenta k x and k y in Eq. (4.1). The main 
physical consequence of such momentum anisotropy is the effective decompactification of 
the Fermi surface, which allows one to focus on a theory with two Fermi patches. Also 
observe that under (4.1) the (d x <p) 2 part of the boson tree level action is irrelevant, which 
justifies omitting this term in eqs. (2.4), (2.7). 

Apart from the fermion time derivative term and the relevant mass perturbation (r — > 
s~ 2 r), all the terms in the Lagrangian (2.7) are marginal. Higher order perturbations to 
(2.7), consistent with the Z 2 symmetry of the order parameter, such as a 4 term, are 
irrelevant. 

We would like to note that for the case of the Ising-nematic (or g-wave) transition the 
low-energy action (2.7) does not possess a <p — > —<P symmetry. This is due to the fact that 
the direction of bosonic momentum q is transformed under ir/2 rotations (reflections) and 
hence the physics is controlled by a different pair of patches of the Fermi surface. Hence, 
in principle, it is possible that in the kinematic regime of interest a 3 term is generated by 
the renormalization group process. Such a term would be marginal under the scaling (4.1). 
A linear term in can also be generated by the effective theory. However, the one-point 
function has momentum q = and, hence, does not belong to any particular kinematic 
regime. In practice, we can demand that the expectation value of <fi is zero in the disordered 
phase by tuning the coefficient of the 0-linear term. In any case, as we will show below, 
there exists a Ward identity, which guarantees that if these terms are initially zero, they 
are not generated by the RG of the low-energy theory (2.7). Note that for the case of the 
spin-liquid or Ising ferromagnet transitions, the low energy theory (2.7) respects the time 
reversal symmetry which maps Fermi patches at ko and — k$ into each other and, hence, 
terms odd in are prohibited. 

An important observation is that the theory (2.7) lacks an expansion parameter. To see 
this, note that due to the rescaling performed in section II, the engineering dimensions, 
[k x ] = [k y ] 2 , but the dimension of to is kept independent. Then, the coupling constant e 2 
has the dimensions [A^] 3 /[u;]. Therefore, e 2 is a dimensionful quantity and cannot be used as 
an expansion parameter. Moreover, e 2 is actually the only parameter in the theory relating 
frequencies and momenta. Hence, its flow under RG is equivalent to an appearance of a 
non-trivial dynamical critical exponent. 

Note that up to this point we have dropped an allowed relevant fermion chemical potential 
term, 

AL = -5^ S (4.2) 
This term can be absorbed into the definition of the momentum fc about which the theory is 
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expanded and, thus, is redundant (note, the scaling dimension [5] = [k x ] = 2). Nevertheless, 
it is convenient to leave this term in the Lagrangian for renormalization group purposes. We 
assume that when the theory is tuned to the criticality r = r c and the coefficient 5 is set to 
5 = 5 C , the Fermi surface passes through the points ko, —ko. 

We now discuss the renormalization of our theory. The Lagrangian contains four marginal 
operators, which each requires a renormalization constant. However, as we will argue below, 
emergent low-energy symmetries of the theory (2.7) imply certain relations between these 
constants. Moreover, the two relevant operators, have the same bare dimension, [r] = [5] = 2. 
Thus, we need to consider possible mixing between these operators. 



A. Rotational Symmetry 

Observe that the initial shape of the Fermi surface does not enter the low-energy theory 
(2.7). In fact, we could have started with a circular Fermi surface with kp = mvF- This is 
reflected by the fact that Eq. (2.7) has an emergent continuous "rotational symmetry", 

</>(x,y) -» (f>(x,y + 0x), ip s {x, y) ->■ e~ is ^ y+ ^ x) ^ s (x, y + Ox) (4.3) 
Equivalently in momentum space, 

q x -9q y -s—,q y + s-j (4.4) 

Note that the rotation angle 9 becomes non-compact and the rotation group becomes M 
instead of U(l). This is a consequence of the effective decompactification of the Fermi 
surface. Moreover, due to the anisotropic scaling 9 is now dimensionful [9] = [k y ]. In 
fact, the situation is analogous to the transformation of the Lorentz symmetry to Galilean 
invariance in the non-relativistic limit u <C c\q\. Here the role of u is played by q x and the 
role of \q\ by q y . 

The symmetry (4.4) implies the following form of the bosonic and fermionic Green's 
functions (we suppress the frequency dependence): 

D(q x ,q y ) = D(q y ) (4.5) 
G s {q x ,q y ) = G(sq x + q*). (4.6) 

In particular, the form of the fermionic Green's function implies that the terms ipl(—isd x )ili s 
and i/j s (—dy)ip s in the Lagrangian (2.7) must renormalize in the same way. Physically, this 
means that the curvature radius of the Fermi surface K does not flow under RG (i.e. K has 
a limit as the cutoff A — > 0). 

The identities (4.5,4.6) ensure that the Green's functions at a given physical momentum 
remain invariant under small changes in the choice of the points ±k on the Fermi surface 
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FIG. 3. The momentum of the fermion at point P can be measured with respect to either the 
co-ordinate system at kg, or that at k\. 

about which the field theory is denned. Let us demonstrate this explicitly using Fig. 3. We 
set the co-ordinate system so that ko = (0, 0), and measure the momentum of a fermion at 
the point P to be (q x , q y ). Now let us shift to the field theory defined at the Fermi surface 
point ki = (k x , Ky). As this point has to be on the Fermi surface, we have k x + = 0. We 
denote the co-ordinates of the point P in the new co-ordinate system by (q' x , q' y ). These are 
obtained from the old co-ordinates by a shift in origin followed by a rotation by an angle 8, 
where tan# = 2n y ; this yields 

q x = q x K X + lKy(C[y Ky) 

q' y = q y -K y , (4.7) 

where we only keep terms to the needed accuracy of 0(x,y 2 ). It can now be verified that 
Qx + l'y = Qx + Qy, and so by Eq. (4.6) the fermion Green's function remains invariant under 
the change in the Fermi surface reference point. Also, by choosing n y = q y we can set q' y = 0, 
and then q x + q 2 is identified as the invariant measuring the distance between P and the 
closest point on the Fermi surface. For the boson Green's function, there is no shift in origin 
of the co-ordinates, and the corresponding transformation is q x = q x + 2n y q y , q' y = q y , and 
this remains invariant under Eq. (4.5). 

These invariances are essential in ensuring the consistency of our description of each pair 
of time-reversed Fermi surface points by a separate 2+1 dimensional field theory. Note that 
such a consistency requirement would not have arisen if we had used a 1+1 dimensional field 
theory at each Fermi surface point, 27 > 28 > 41 ~ 45 because then every fermion momentum would 
appear only in the theory denned at the closest point on the Fermi surface. In our case, we 
are free to use the 2+1 dimensional theory at this closest point, or at any of the neighboring 
points. 
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Before concluding this section, we would like to point out that in the case of the Ising- 
nematic transition, the "rotational symmetry" (4.4) is not related in any way to "large" 
rotations by 7r/2, which are actually not implemented in the low-energy theory. 

B. Ward Identities 

We now examine the consequences of Ward identities associated with the global symme- 
tries of Eq. (2.7). Similar consequences were implicit in the analysis of the superconducting 
case in Ref. 36. Here we will present a more formal analysis, which also shows that Eq. (3.20) 
in Ref. 36 holds to all orders in 1/N. 

The low energy theory (2.7) has two continuous global U(l) symmetries. The first of 
these is related to the conservation of particle number, 

U(1) F : tp+ -> e ia i>+, V>- -> e*>- (4.8) 

The conserved current associated with this symmetry is, 

Ur,jx,j y )F = (iv(^W+ + V-il>-),V+il>+ ~ V-i>-, -i{il>\. d y i)+ + d (4.9) 

For the spin-liquid problem, the gauge field couples precisely to the x component of jf- 

The second U(l) symmetry is lattice translation. Indeed, ip + and come from opposite 
points in the Brilloin zone and, hence, transform under general lattice translations as, 

U(1) T ■ ^+ -> e ia ip+, ^_ -> e _ia V- ( 4 - 10 ) 
The conserved current associated with this symmetry is 

(jr,3x,3 y )T = (irj(iplip + - + ip-ip-, -i(il>+ 9 - ^1 d (4.11) 

Observe that the Ising-nematic order parameter couples to the x component of Jt- Note 
that despite the similarity of the spin-liquid and Ising-nematic problems, there is an impor- 
tant difference. In the spin-liquid case, the gauge field couples to the fermion current on all 
energy scales. In the case of the Ising-nematic transition, the order parameter couples to a 
conserved current only at low energies. 

We note in passing that for an Ising ferromagnet transition, the current to which the 
order parameter couples is related to the symmetry, 

U(l)! : V +t -+ V-t^e- fa ^_, rf, u e~^ +h ^ -»■ e ia ^ (4.12) 

In fact, this is not a symmetry of the underlying theory, but only of the low-energy La- 
grangian (2.4). The symmetry is broken by four- Fermi interactions, which are however 
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irrelevant under (4.1). 



Current conservation implies that the insertion of d T j T + d x j x + d y j y into any correlation 
function is zero, up to contact terms (we have dropped the current subscript; the current, 
which couples to the order parameter is implicitely assumed). We note that the temporal 
component of the currents (4.9), (4.11) has a coefficient r\ in front and, therefore, can be set 
to zero in the kinematic regime of interest. We, thus, have d x j x + d y j y ~ 0. Defining the 
one-particle irreducible polarization function, 

HiiGz) = J drd 2 xe^ T -^{ 3l {x) 3j {Q)) 1PI (4.13) 

we have 

q x U xx {q) + q y Il yx (q) = (4.14) 
We note that Ii xx (q) = U xx (q T ,q y ) is precisely the irreducible boson self-energy. Hence, 

q y 

Power counting indicates that U xx has the following UV structure 

IW<? T , q y ) U = K x + K 2 r + K 3 q 2 y (4.15) 

where K\ ~ A 2 , K2,K^ ~ log A and A is the UV cut-off with dimensions of q y . For 
U yx (q T ,q x ,q y ) to have an analytic UV behaviour (as again expected from power counting), 
we must have 

K 1 = K 2 = 

Thus, the coefficient of the mass operator <f) 2 requires no renormalization (i.e. the metric 
factor Z r has a limit as A — > 0). 

An interesting question is whether the polarization function U xx actually vanishes for 
q y — > as suggested by Eq. (4.14). However, for finite q T we already know from one-loop 
calculations that such a limit does not exist within the scaling regime, as 

n xx (q T , q y )iioop Q>"i T; n yx (q T , q Xj q y ^n 00 p c^- 



\q y \ q y \q y \ 

However, one might hope that the limits \\m qy ^ Q lim^^o n xx (q T , q y ), U xy (q T , q x , q y ) do exist. 
In this case, we would conclude, 

lim lim Il xx (q T , q y ) = (4-16) 
which would be a stronger statement than the non-renormalization of the mass term. Oth- 
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erwise, if the limit above exists only for U xx by not H xy then, 

lim lim U xx (q T ,q y ) = c r r (4.17) 

with c r - some universal constant. We have explicitly checked that to three loop order c r = 
and the strong form of the non-renormalization identity Eq. (4.16) holds. 

One can generalize the discussion above to higher order correlation functions of the order 
parameter. Ward-identities imply that the effective potential for the (ft field is not renormal- 
ized from its tree-level form, 

V((ft) = T -<P 2 (4.18) 

This property is also shared by the theory of the nematic transition in a <i-wave superconductor. 35 ' 36 
In particular, no (ft 3 term is induced in the Lagrangian by the renormalization group process 
if this term is originally zero. (Note that if a (ft 3 term is initially present, correlation func- 
tions of currents no longer coincide with the correlation functions of the order parameter, 
and the Ward identities do not constrain the renormalization properties of the theory). The 
effective potential (4.18) becomes unstable for r < 0. Thus, we expect that in the ordered 
phase the theory is controlled by dangerously irrelevant operators, such as (ft 4 . 
Finally, one can derive a Ward identity for the fermion boson vertex, 

q x T x (q,p,p + q) + q y T y (q,p,p + q) = G^(p + q) - G~ 1 (p) (4.19) 

with 

T i (q,p,p + q) = J cfe r d 2 £^ r dV~^' T+ ^ (4.20) 

G(p) = [ drd 2 xe iqTT -^ s {ift{x)ift\0)) (4.21) 



T x is precisely the irreducible fermion-boson vertex. Power counting gives UV structure of 
T x and G~ l as, 



r x {q,p,p + q)=C 1 (4.22) 
G-\p)=C 2 + C 3 ( Px +p 2 y ) (4.23) 

Thus, for the UV behaviour of T y to be analytic in external momenta, C\ — C3. Therefore, 
the vertex and the fermion self-energy renormalize in the same way. Hence, the boson field 
requires no field-strength renormalization (i.e. the metric factor has a limit as A — > 0). 

Before concluding this section, we would like to note that perturbation theory based on 
self-consistent propagators (3.4), (3.6) actually does not respect the Ward identities. This is 
due to the fact that these one-loop propagators include the fermion self-energy correction, 
but not the vertex correction. However, since the fermion self-energy is only frequency 
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dependent, Ward identities involving currents at zero external frequency are still respected. 



C. RG equations 

From the discussion above, we conclude that at criticality, our theory needs only two 
renormalizations: a rescaling of the field strength of the fermion field i/j and a renormalization 
of e 2 , 

V = 4 /2 </v, e 2 = Z e e 2 r (4.24) 

Here the subscript r denotes renormalized quantities and we define renormalized irreducible 
correlation functions of boson and rif fermion fields as, 

T n b ,n f = z n f /2 rnbtnf ^ ^ 

Both Z^ and Z e are functions of A//x where \l is a renormalization scale (which we choose 
to have dimensions of q y ) and of the number of fermion flavours N. As e 2 is dimensionful, 
Z^p and Z e cannot depend on it. We introduce the anomalous dimensions, 

b = A-^\ogZ e (4.26) 

7to = -A^log^ (4.27) 

The constants t]^ and b are expected to be pure universal numbers, independent of A//i. 

Away from criticality, we recall that by the Ward identity, the coupling r does not renor- 
malize. On the other hand, the coupling 5 can pick up a renormalization linear in r, 

S = S c + S r + Z rS e 2 r r (4.28) 

with Z T & again a function of A//i only. In what follows, we denote 5 — 5 C as S for brevity. 
Note that there is no renormalization constant in front of S r since a finite change in 5 only 
shifts the value of k x in correlation functions: 

r n "' n f({p},S + a) = r n »' n f({p-sax},5) (4.29) 

where s = ±1 for momenta of fermions ip± and s = for boson momenta. We let, 

a = Z^A-^Zrs (4.30) 
Now, differentiating Eq. (4.25) we obtain the renormalization group equations 

( A M + be ^ + ae2r §s - fy*) rW {{Pyh {Pxh M ' r ' 5 ' e2 ' A) = (431) 
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It is convenient to get rid of the derivative with respect to 5 in Eq. (4.31). To do so, let 
the location of the Fermi-surface of fermion at finite 5 and r be given by k x + k 2 = 
Ak(r,5,e 2 ,A). Then, Ak is clearly a physical quantity and must satisfy, 

a|t + + Ak(r, 5, e 2 , A) = 0. (4.32) 



OA de 2 d5, 
We will solve this equation shortly. However, first note that 

dAk 



as L < 433) 

Now, it is convenient to expand momenta around the physical Fermi- surface, defining, 

fw({p},r,5,e 2 ,A) = r nb ' nf ({p + sAk(r, 5, e 2 , A)x}, r, 5, e 2 , A) (4.34) 
The resulting T is independent of 5 and by Eqs. (4.29), (4.31), (4.32), (4.33) satisfies, 

( A Jv +be2 ^ ~ if**) r n6 ^(K},K},{^},r,e 2 ,A) = (4.35) 
By dimensional analysis, 

and solving the RG equation, we obtain 

f n ^(s{p y }, s 2 {p x }, s 3 ~ b {u}, s 2 ~ b ?) = ^H^/^/^ {p x .}, {a;}, f) 

(4.37) 

Hence, the critical theory is invariant under, 

p y ->■ pa ->■ s 2 p x , w -> s^u; (4.38) 

with 

z = 3-b, (4.39) 

where z is the dynamic critical exponent. Note that we have defined z with reference to 
length scales associated with directions tangent to the Fermi surface (y); as indicated in 
(4.38), length scales orthogonal to the Fermi surface scale as the square of length scales 
tangent to the Fermi surface. Moreover, if we define £ as the correlation length along the y 
direction then upon approaching the critical point, £ ~ r~ u , with 

v = — !— -. (4.40) 

Z 1 
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Note that by combining Eqs. (4.24,4.26,4.39) we can write down the RG equation for the 
coupling e: 

a mL, = -(*- 3 ^ (441) 

This shows that the renormalization of the coupling e is directly related to the dynamic 
critical exponent, as we had claimed earlier. 

Now, let us consider a few explicit examples of correlation functions. For the bosonic 
two-point function we have, 



D-\q y ,u) = rg (q y {re 2 A z ~ 3 )-^ , u(r 2 e 2 r 3 )^ 



(4.42) 



Note that, 



lim lim D l (q v ,uj) = rg(0,0) 

q y ->0 w^O 



i.e. the Ising-nematic susceptibility satisfies \ ~ r 7 with the exponent 



(4.43) 



7 = 1. 



(4.44) 



We may also write more succinctly, 

D-\q y , u>) oc r {z - 1} 9(q y C, coe 2 A z - 3 C) (4-45) 



So far, we have been concentrating on a fixed direction of bosonic momentum q. Now let 
us study the dependence of the result on q. Using Eq. (2.6) 

D-\q,u) = ZfK-^rog (igK^VA^ro)"^, u;(Z r - z e 2 A^ 3 r z )-^i) (4.46) 

where for brevity tq is taken to denote the deviation from the critical point. We concentrate 
on the static limit u — 0. In a Fermi liquid, the susceptibility must have a continuous limit 
as q — > 0. Therefore, we conclude that the combination Z^KZ r must be independent of the 
direction q. This is quite plausible, as neither of the constants run under RG. 
Now let us look at the behaviour of susceptibility at the critical point, 

In particular, the static susceptibiltiy satisfies, 

D-\&U) ~ aftM*- 1 (4.48) 

In the context of the spin-liquid problem, many studies 50 ' 52-56 examined the structure 
of the higher loop corrections to the susceptibility. In particular, Kim et al. 50 examined 
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two-loop corrections to Im Z} _1 (<f, a;) for real frequencies \u\ <C \q\, and found no corrections 
to the leading answer ~ w/\q y \ in Eq. (3.4); Fermi liquid arguments were made 50,52-54,56 that 
this functional form held at higher orders. However, this result by itself does not fix the 
value of z; indeed, ImD _1 (g, u) ~ u/\q y \ is consistent with the scaling form (4.47) for any 
z. These studies also implicitly assumed a Fermi liquid picture with D -1 (<f, u = 0) ~ q 2 , 
and this does imply z — 3. We will examine _D _1 (<f, u = 0) up to 3 loops in Section VA, 
and find no correction to z = 3. 



Proceeding to the fermion Green's function, 

Gj\k,cu) = A 2 (^j L (fc(re 2 A- 3 )-^ lU (rVA- 3 )-^) (4.49) 



with k = sk x + ky - the distance to the Fermi surface. More compactly, 

G-\k, u) oc £-( 2 "^)L (kg*, ue 2 A z - 3 £ z ) (4.50) 

A crucial property of the theory that is manifested by the above expression is that the 
"fermionic correlation length" scales as the square of the "bosonic correlation length" . 



For u <C £ z , k <C £ 2 we expect the fermion Green's function to assume a Fermi-liquid 
form, 

G(k,u)= . j - (4.51) 
— too + VFk 

By matching to the scaling form, 



£-( z - 2 ), Z ~ £"(^- 2 ) (4.52) 



Notice that both the Fermi velocity Up and the residue Z tend to zero as we approach the 
critical point, albeit with different power laws. Finally, at the quantum critical point, 

G-\k,u) = A^k'-^P (j^^j , (4-53) 

where we reiterate that k = sk x + k 2 is the distance to the Fermi surface. In particular, the 
self-energy on the Fermi surface scales as, 

G~\0,u;) ~ J 2 -^ )/z (4.54) 

and the static self energy, 

G~\k, 0) ~ fc 1 "^/ 2 (4.55) 
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Moreover, from Eq. (4.53) we can obtain the tunneling density of states, 



N(u) = 




(4.56) 



where 

A(k,w) = —lmG(k,iu -> u + iO + ) (4.57) 

7T 

The k integral in Eq. (4.56) factorizes into integrals over components along and perpendic- 
ular to the Fermi surface. The former gives a factor proportional to the perimeter of the 
Fermi surface, while the later yields the frequency dependence, 

N(u) ~ d** 1 * (4.58) 

We remind the reader that the expression in Eq. (4.58) corresponds to the physically ob- 
servable electron tunneling density of states only in the case of a nematic transition, as for 
the spin/charge-liquid problem, the physical electron operator is a product of if) and a boson 
operator. 

Related scaling forms for the fermion Green's function were discussed on a phenomeno- 
logical basis by Senthil. 58 However his definition of z differs from ours. We define it using 
the fermion momentum parallel to the Fermi surface, because this is the natural momentum 
scale appearing also in the boson correlations. He defines it by the fermion momentum 
orthogonal to the Fermi surface, which scales as the square of the parallel momentum. 

Finally, let us discuss the shift of the Fermi surface Ak. Using Eq. (4.33) in the RG 
equation (4.32), we obtain, 

Ak = -?—re 2 + C k (re 2 A z ~ 3 ) 2u + $ (4-59) 
z — 3 

Thus, the shift of the Fermi surface upon deviation from the critical point receives two contri- 
butions: one analytic in r and the other non-analytic. Reexpressing the second contribution 
in terms of the correlation length, 

Ak = -re 2 + C k C 2 + $ (4.60) 

Z O 

where the coefficient C k is expected to be universal. We would like to point out that the 
case z = 3 has to be treated separately. In this situation one obtains, 

2 2 
OITP TP a 

Ak = bg JL = - Ck C 2 log(A£) + 5 (4.61) 

with C k again universal. 

The value of the Fermi surface shift Ak can be used to compute the compressibility, I 2 , 
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where /x is the physical chemical potential. Indeed, by Luttinger's theorem the change in 
density can be obtained as, 

n r 

Sn = — — / dsAk(9) (4.62) 



(2tt) 2 j 

where the integral is over the circumference of the Fermi-surface. The main question is how 
does the chemical potential enter our low-energy theory. If /x only couples to the operator 
ip^ip, renormalizing the value of 5, then from Eqs. (4.59), (4.61) we would conclude that 
the compressibility tends to a constant and has no interesting corrections near the quantum 
critical point. On the other hand, if the coupling r has a non-trivial /x dependence, then we 
would conclude, 

— = k + mC- 3 , z ± 3 (4.63) 

— = Ko + Kik>gAf, z = 3 (4.64) 

Note that for z > 3 the above forms imply that the compressibility diverges as we approach 
the critical point. 



V. ANOMALOUS EXPONENTS TO THREE LOOPS 



In this section, we evaluate the exponents z and r]^ to three loop order. We find that the 
exponent rj^ is non-zero at this order. The value of rj^ is not suppressed in the large- N limit. 
On the other hand, the dynamical critical exponent z remains unrenormalized from its RPA 
value z = 3 to this order. Moreover, in the large- N limit, the boson self energy acquires a 
finite correction of order N 3 ^ 2 , which is larger than the bare value (order N). Finally, we 
find that the constant a in Eq. (4.30) associated with the shift of the Fermi-surface away 
from criticality is non-zero at three loop order. We note that the N 3 ^ 2 correction to the 
boson self-energy and the non-zero rj^ are only present for the Ising-nematic and spin-liquid 
universality classes, and are absent for the Ising ferromagnet transition. 



A. Dynamical critical exponent 

Let us first address the question of renormalization of e 2 . At two- loops the only correction 
to the static boson-self energy U(q T = 0,q), which is not already taken into account by the 
solution to self-consistent Eliashberg equations is given in Fig. 4. However, this diagram 
vanishes when the external frequency is equal to zero. Indeed, as pointed out in Ref. 40, any 
diagram with fermions from a single patch, in which the fermion propagators involve a sum 
of two or less internal momenta, vanishes in the static limit (one picks the internal frequency 
with the largest absolute value and integrates over the corresponding x component of the 
momentum. All poles will be in the same half-plane). Actually, a calculation presented in 
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Appendix B shows that the diagram in Fig. 4 vanishes for any external frequencies and 
momenta. 




FIG. 5. Three loop corrections to the boson self-energy with one fermion loop. 

The three loop corrections to n(g) are shown in Figs. 5 and 6. By the argument described 
above, all of these diagrams vanish when the external frequency is zero if all the fermions 
are from the same patch. Hence, the only non-zero corrections to Il(g T = 0, q) come from 
the Aslamazov-Larkin type diagrams, Fig. 6, 
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(b) 

FIG. 6. Aslamazov-Larkin type three loop contributions to the boson self-energy. 



S 3 U(q) 



dl T d 2 l 



r 3 (g, l,-{l + q))T 3 (-q, -/, / + q)D(l)D(l + q) 



(5-1] 



Here T 3 is the fermion-induced cubic boson vertex, which receives contribution from the two 
fermion patches, 

r 3 = r 3 + r 3 (5.2) 
r1(/i, h, h) = N\ 3 s (f s (h, z 2 , i 3 ) + f s (i 2 , h, k)) 

dp T d 2 p 



fs(h, hi h 



-G s (p)G s (p-h)G 3 (p + h) 



(5.3) 
(5.4) 



The diagrams where the fermions in the two loops come from the same patch give a vanishing 
contribution to H(q T = 0, q). Thus, to three loops, 



5 3 U{q T = 0,q) 



dLd 2 l 



(27T 



3 T 3 + (q, /, -(/ + q))fi(-q, -I, I + q)D(l)D(l + q) + (q 



-q) 



-A 3 A 3 A^ 2 



dl T d 2 l 



f + ( qj l,-(l + q))(f_(-q,-l,l + q) 



+ f4-q, I + q, -l))D{l)D{l + q) + (q -> -q). 



(5.5) 



The two terms in brackets in the equation above originate respectively from diagrams in 
Figs. 6 a) and b). Converting these diagrams into the double line representation of Ref. 40, 
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(a) 




T 




(b) 



FIG. 7. Double line representation of Ref. 40 applied to the Aslamazov-Larkin diagrams in Fig. 6. 
The fermions in the two loops are assumed to come from opposite patches. We have reversed the 
directions of the fermion propagators from the second patch, and the dotted arrows indicate the 
true directions of the fermion momenta. 



we obtain Figs. 7 a) and b). [We remark that the genus expansion of Ref. 40 was developed 
for a theory with only a single Fermi-surface patch. The extension to the present case of 
a pair of time reversed patches is simple: a reversal of the direction of loops with fermions 
from the second patch reduces the problem to that with one patch only. The diagrams in 
Fig. 7 have their lines reversed precisely in this way. The additional dotted arrow besides 
each propagator indicates the true direction of fermion momentum.] In this representation, 
the graph a) contains a loop while the graph b) does not. As a result, in the genus expansion 
of Ref. 40, the diagram in Fig. 6 a) is enhanced to O(N), while the diagram in Fig. 6 b) 
is of 0(1). However, we will see that the diagrams are actually individually ultra-violet 
divergent, as a result the counting of Ref. 40 is inapplicable here. It turns out that the sum 
of the diagrams is UV finite and of 0(N 3/2 ). 

We give details of the evaluation of Eq. (5.5) in Appendix B, where we find 



U(q T = 0,q) = C\ + \ 




(5.6) 



In the large- N limit, the coefficient C is given by, 



C « -0.09601N 3/2 , N — > oo 



(5.7) 
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while for the physical value N = 2, 

C « -0.04455, N = 2 (5.8) 

The iV 3 / 2 behaviour in Eq. (5.7) indicates a breakdown of the genus expansion of Ref. 40. 
Moreover, since this correction is parametrically larger than the tree level value, the existence 
of the large- iV limit of the theory is cast into doubt. In particular, it is not clear if there 
are higher loop graphs with even stronger divergences in the large-iV limit. Moreover, we 
expect contributions to the bosonic self-energy analytic in q y to be generated from kinematic 
regimes involving the whole Fermi-surface and not just the two Fermi patches. Such analytic 
contributions might also exhibit anomalous scaling with N. 

Note that there is no logarithmic dependence on A/fi in Eq. (5.6), and so we have z = 3 
at this order. For the physical value of N = 2, the finite three-loop correction turns out to 
be rather small numerically. 



B. Fermion anomalous dimension 

The Feynman diagrams for the fermion self-energy up to three loop order are shown in 
Figs. 8, 9 and 10. By reasons explained in the previous section, the diagrams in Figs. 8 and 

9 vanish when the external frequency is zero and, hence, do not contribute to the fermion 
anomalous dimension. 

/ / \ \ 

> L » i „ 

P p-h p-h-h p-k P 

FIG. 8. Fermion self-energy at two loops. 

Thus, the only fermion self-energy diagrams that can give UV divergences are shown in 
Fig. 10. Actually, the diagram in Fig. 10 a) is zero since the polarization correction in Fig. 
4 vanishes. Thus, we only need to consider the two diagrams in Fig. 10 b) and c). For these 
graphs to be UV divergent, the fermions running in the loop and the external fermions must 
come from different patches. The diagram in Fig. 10 b) contains two loops in the double 
line representation (Fig. 11 a)) and is expected to be of order 1/N, while the one in Fig. 

10 c) has no loops in the double line representation (Fig. 11 b)) and, hence, is expected to 
scale as 1/N 2 . 
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/ / \ \ \ 

, / , / , r + ! 2 \ l \ \ , \ , 

P p + h P + h p + h-hp + h P 



> I > ' > I > \ > ^ > \ > 



'V N 

/ h-h \ 

„ i t / , p\ l n\ i > \ > 

P p + h p+h P + h-hp-h P ; 



, / , / , \ p ;' 2 / , \ , \ , 

P p- hp- h-h p-h-hp-h P 



I h 11 

/ / /\ \ \ 
/ / / \ \ \ 

P p- hp -h-h p-hp+h-h P 



FIG. 9. Three loop fermion self-energy diagrams with no fermion loops 



A calculation presented in Appendix B gives the UV divergent contribution, 

5 3fe S + (c = 0,p) = X + X-J b ( Px + pl)log ( ^ + A ^ |1/2 ) , (5.9) 
5 3c £+(u; = 0,p) = 5 3c Z + (cu = 0,p = 0) + A + A_ J c (p x + p\) log ( \ 

V \\Px+Py\ 1/2 

The constant Jb is independent of N and given numerically by, 



(5.10) 



J b ^ 0.1062 (5.11) 
On the other hand, the constant J c is ^-dependent. For N = 2 we obtain, 

J c « -0.03795, iV = 2 (5.12) 
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(a) 



i 



k k 




V p-h P-h P P p + h p+h P 

(b) (c) 



FIG. 10. Three loop fermion self-energy diagrams with one fermion loop 




(b) 

FIG. 11. Double line representation of fermion self-energy diagrams in Figs. 10b, c, as in Fig. 7. 
The external fermions and the fermions inside the loop are assumed to come from opposite patches. 

while in the large- N limit, 

Jc "i^ l0g3iV ' iY ~ >0 ° (5 - 13) 

Notice that there is no 1/N suppression in Eq. (5.9). A way to interpret this, is that the 
diagram is really of order 1/N (as the genus expansion predicts), however, it is a function 
of N{p x +Py)- Indeed, recall that the genus expansion assumes N{p x + Py) ~ 1. However, 
the UV divergent piece of the diagram cannot depend on the magnitude of p x + p 2 y and is 
valid for any external momentum or frequency. On the other hand, the infrared scale under 



27 



the log is expected to become u 1 ^ 3 once oj ^> N 3 / 2 \p x + p 2 \ 3 ^ 2 - Also observe that up to a 
logarithmic enhancement, the non-planar diagram 10 c) (11 b)) is of order 1/N 2 , as expected 
from the genus expansion. 

Note that the UV divegence in Eqs. (5.9), (5.10) is logarithmic, as expected from power 
counting, and comes from a region where both internal momenta and frequencies diverge 
in accordance with the scaling (4.1). This is unlike the anomalous linear divergences of 
the Aslamazov-Larkin diagrams that occur when the internal momenta q y are of order of 
external momenta, while internal q x , q T diverge. 

Thus, to three loop order, 

5 3 Z + (tu = 0,p) = 5 3 Z + (tu = 0,p = 0) + X + X-J( Px +p 2 y ) log ( - ^ 2|1/2 ) (5.14) 



, 0.06824 N = 2 . 1 . 

J = J b + J c ^ { (5.15) 
* 0.10619 N = oo ( J 

Although the self-energy correction (5.14) is not parameterically suppressed compared to 
the bare value even when A" = oo, it appears to be suppressed numerically. Thus, we may 
estimate, 

Z^ — l — A + A_JlogA//i 

^ = \ + \_J = ±0.06824 (5.16) 

where the upper sign refers to the Ising-nematic transition and the lower sign to the spin- 
liquid and we have used the value of J at N = 2. 



C. Fermi surface shift 



We now evaluate the coefficient a, Eq. (4.30), associated with the renormalization of 
chemical potential 5 away from criticality. This coefficient can be obtained from the insertion 
of the <j) 2 operator into the two-point fermion Green's function at criticality. By setting 
all external frequencies to zero, we find that at three loop order the only UV divergent 
contribution can originate from the diagrams in Figs. 10 b) c) with the 2 operator inserted 
into the boson propagators. The details of the calculation are presented in appendix B. We 
find, 

^ = Jre 2 log A y (5.17) 
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with 



J r = 0.00208, N = 2 




(5.18) 



Absorbing this divergence into the chemical potential, 



Z rS = J r log Ay/ fl 



(5.19) 



and 



(5.20) 



Thus, the <fi 2 operator mixes with the i/j^ip operator. If the dynamical critical exponent 
z = 3, this leads to a logarithmic divergence of the compressibility, Eq. (4.64). Note that 
the magnitude of the mixing a is suppressed in the large N limit and is also numerically 
small for N = 2. 

VI. CONCLUSIONS 

This paper has presented the scaling properties of the field theory in Eq. (2.7) which 
describes a number of problems involving the breakdown of Landau Fermi liquid theory at 
all points on a two dimensional Fermi surface. The main motivation was provided by the 
quantum phase transition caused by the onset of Ising-nematic order, which reduces the 
point-group symmetry from square to rectangular. However our theory also directly applies 
or can be generalized to breaking of other point-group and/or time-reversal symmetries, 
and these were described in Section II. One of these cases is the "circulating current" 
order parameter of Simon and Varma 33 ' 47,48 . Apart from applications to quantum critical 
points, our theory also described non- Fermi liquid phases associated with spin liquids 57 ' 58 or 
algebraic charge liquids 59-61 , which have Fermi surfaces coupled to U(l) gauge fields. 

Our critical theory was formulated in terms of a time-reversed pair of patches on the Fermi 
surface, centered at the wavevectors ±k (see Fig. 1). The value of k was determined by 
requiring that the tangent to the Fermi surface at ko be parallel to the wavevector q carried 
by the order parameter insertion in the correlation function being computed. However, in 
general, there is nothing special about the point ko, and neighboring points on the Fermi 
surface should behave in a similar manner. This key feature was implemented in our theory 
by the rotational symmetry discussed in Section IV A, and the identities (4.5,4.6), which 
show that the Green's function remains invariant as we move along the Fermi surface. 

We emphasize that although we have critical theories associated with every pair of points 
on the Fermi surface, the Lagrangian (2.7) and all the fields are 2+1 dimensional i.e. <j) and 
ip a are integrated over arbitrary functions of x, y, and r. Thus, as we noted earlier, our 



29 



approach and results differ from studies using a 'tomographic' representations of the Fermi 
surface, in which every point on the Fermi surface is described by a 1+1 dimensional field 
theory. 27 ' 28 ' 41-45 Our 2+1 dimensional representation leads to a redundancy in our description 
of the degrees of freedom, and the identities of Section IV A ensure the consistency of this 
redundant description. 

Our main results include the scaling relations for the order parameter susceptibility in 
Eq. (4.45), and for the fermion Green's function in Eq. (4.50). These are associated with only 
two independent exponents, the dynamic scaling exponent z, and the fermion anomalous 
dimension r\^. The correlation length exponent v was given by exact scaling relation in 
Eq. (4.40), while the susceptibility exponent 7 = 1. For the spin-liquid case, Fermi liquid 
arguments were m ade 49 ' 50,52 ~ 54 ' 56 suggesting that z = 3; we found z = 3 to three loop order in 
Section V, although we did not prove this to all orders, and our scaling theory is compatible 
with a general value of z. Our three loop computation also gave a non-zero value of 77^, 
with opposite signs for the Ising-nematic and spin-liquid cases. In the case of the nematic 
transition, a non-zero positive 77^ implies the suppression of the electron tunneling density 
of states, Eq. (4.58). Another striking effect that we find for the case of a nematic transition 
is the power law divergence of the compressibility for z > 3, which turns into a logarithmic 
divergence if z = 3. 

Our scaling results were expressed in terms of correlators of the fermionic field ip +a 
carrying momentum q as measured from the point k from the Fermi surface, implying from 
(2.3) that the electron c a has momentum ko + q (and similarly for ip- a ). However, note that 
(after appropriate rescaling of momenta, and for a circular Fermi surface) \k\ — &f ~ q x + g 2 . 
Thus the identity (4.6) implies that the scaling function (4.50) for the two-point fermion 
Green's function depends only on \k\ — kp. This is similar to the dependence found in 
other treatments e.g. in the recent critical theories 62 ^ 65 obtained by applying the AdS/CFT 
duality to fermions propagating near a Reissner-Nordstrom black hole. The latter theories, 
in their current classical gravity formulation, find 64 = 0. 

It is also interesting to compare the structure of the critical theory in the AdS / CFT frame- 
work to that found here. We have an infinite set of 2+1 dimensional field theories labeled by 
pairs of momenta on a one-dimensional Fermi surface i.e. a S 1 /^ set of 2+1 dimensional 
field theories. In the low-energy limit, the AdS/CFT approach yields 64 a AdS2 x 1R 2 geom- 
etry: this can be interpreted as an infinite set of chiral 1+1 dimensional theories labeled by 
a M 2 set of two-dimensional momenta k. It is notable, and perhaps significant, that both 
approaches have an emergent dimension not found in the underlying degrees of freedom. 
We began with a 2+1 dimensional Hamiltonian, and ended up with a § 1 /Z 2 set of 2+1 
dimensional field theories. In AdS/CFT, there is the emergent radial direction representing 
energy scale. These emergent dimensions imply redundant descriptions, and require associ- 
ated consistency conditions: we explored such consistency conditions in Section IV A, while 
in AdS / CFT the consistency conditions are Einstein's equations representing the renormal- 
ization group flow under changes of energy scale. It would be interesting to see if fluctuations 
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about the classical gravity theory yield corrections to the A0IS2 x M 2 geometry which clarify 
the connection to our theory. 

In the analysis of the spin-liquid problem, Ref. 40 considered a single patch of the Fermi 
surface, and argued that the 1/N expansion should be organized by the genus of the Feynman 
graph (after the propagators are written in a suitable double line representation, and the 
graph is interpreted as lying on a two-dimensional surface). In our two-patch theory here, 
we have shown that this genus counting is violated. This is the implication of the A 3 / 2 
dependence of the boson self-energy in Eq. (5.6). In fact, at present, it is not clear how 
to take the large- N limit of the theory. On the other hand, for the physical value N = 2, 
we found that the higher loop contributions are numerically small, which suggests that the 
critical exponents are close to the Hertz mean-field values. However, because the loop-wise 
expansion does not possess even a formal expansion parameter, it is not clear if there is a 
systematic way to extract corrections to the mean-field exponents. Thus, our value of the 
fermion anomalous dimension r]^, Eq. (5.16), should be regarded as an estimate only. 
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Appendix A: Decoupling of non-collinear momenta 

In this section we will argue that the fluctuations of the order parameter at non-collinear 
momenta effectively decouple. We focus for simplicity on the case of an Ising-nematic 
transition. We follow the standard Hertz approach, integrating out the fermions to obtain 
an effective action for (f), 




d x 1 ...d x n T n (x 1 ,x 2 , . . .x n )(j)(xi)(j)(x2) ■ ■ .<fi(x n ) 



(Al) 



The n-point effective vertex T n is given by, 



r n (q 1 ,q 2 ,...q n ) 



N 



/"(<?!, <?2, •• • Qn) + permutations of q 1 ,q 2 ,...q. l 



(A2) 



11 



n 



with 



r r\h d 2 k ^ r 

re*, • • • *o = J n p + ^54* 



(A3) 
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where U = Y^j=i Qj- For now we work with "undressed" propagators, 



G(u, k) = — l -— — (A4) 

with k - the distance to the Fermi surface and vf{0) - the local Fermi velocity. As is 
well-known, for oj <C Vp\q\ and \q\ -C kp the two-point vertex has a Landau-damped form, 



T 2 (u,q) = N 



II -Q, 

7(9) -pr + -o + r 
91 e 2 



(A5) 



where the coefficient of the non-analytic term 7(g) = Kd 2 /(2irv F ) with the Fermi-surface 
curvature radius K, Fermi- velocity Vp and form-factor d evaluated at the point on the Fermi 
surface to which q is tangent. On the other hand, the coefficients of the analytic terms r 
and 1/e 2 come from the entire Fermi- surface. 



If we truncate the series (Al) at the quadratic order 

N f dud 2 q 



S 2 = -^ 



II -Q 

1 91 e 



\m«>)\ 2 (A6) 



2 J (2tt) 3 

then at the critical point r = the action (A6) is invariant under the scale transformation, 

<f>(x, t) -> s 3/2 0(sf, s 3 r) (A7) 

Note that here, in contrast to Eq. (4.1), all components of q are scaled in the same way as we 
are not studying the effects of fluctuations with collinear wave- vectors. We can regard the 
terms in Eq. (Al) with n > 2 as perturbations to the Hertz action (A6). Hertz noted that if 
the effective vertices T n possess a regular expansion in frequencies and momenta, such that 
the corresponding operators can be represented as polynomials in the order parameter </> and 
its derivatives, then the perturbations with n > 2 are irrelevant due to the large effective 
dimensionality, -D e // = d + z — 5, with d = 2 - spatial dimension and z = 3 - the dynamical 
critical exponent. Indeed, the perturbation J d 2 xdr(j) n (x) scales as s 3 "-/ 2-5 under (A7) (in 
the special case n — 3, the operator <fr 3 is actually prohibited by the 90° lattice rotation 
symmetry. The lowest dimension local operators with three powers of <ft that are allowed 
by symmetry are 4>((d x (j)) 2 — (d y (f)) 2 ) in the d x 2_ y 2 case and (pd x (pd y 4> in the d xy case, which 
scale as s 3//2 ). 

However, due to the presence of low-energy excitations on the Fermi-surface there is no 
reason to expect that the effective vertices T n would possess a regular expansion in momenta. 
Indeed, we have already seen that the two-point vertex has the non-analytic Landau-damped 
form (A5). As we now show, similar non-analyticities occur in the higher order vertices. 

Let us estimate the vertices (A2) when the external frequencies and momenta obey the 
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Hertz scaling (A7), w~ |g| 3 , <f — >■ 0. In this regime, 



f n (q 1 ,q 2 ,. . .q n ) 



dk T dkd9 
(2tt) 3 



dki 



dd 



n-1 



d(er I 



1 



i=0 



-i(k r + l iT ) + v F (9)(k + v F {9) ■ h 



(A8) 



Let us perform the integral over k. Observe that if |fc r | > Q, with Q = maxj \ U T \ then the 
integral vanishes as all the poles of the integrand are in the same half-plane. Thus, the range 
of the internal frequency is limited by the external ones. With this in mind, 



f n (qi,q 2 , ■■■qn) 



\k T \<n 



dk T 

n-1 



dd 

2^ 



dkt 



dd 



d(9) n 

Me) 



n-1 



3=0 i=0,i¥=j 



4(l iT - l jT ) + V F (9) ■ (U - lj) 



(A9) 



where we have used the symbol for the step function, to avoid confusion with the angular 
variable 9. Now, since q T ~ |<?| 3 /(7e 2 ) v F \q\, for general 9 we can ignore the frequency 
dependence in the denominator of Eq. (A9). Then the angular integration yields a factor of 
0(1) and the integral over k T yields a factor of external frequency, so that 



r n (gi,g 2 , ■■■q n ) 



(A10) 



Note that the momentum dependence in Eq. (A10) is far from analytic. Also, note that for 
n = 2 the result is consistent with the standard Landau damping. 

The only possible caveat to the estimate (A10) is associated with regions of angular 
integration where v F (9) ■ (U — lj) — > 0, i.e some combination of external momenta becomes 
tangent to the Fermi surface. Then the angular integration acquires poles just off the real 
axis, with the imaginary parts of the poles provided by the frequency dependence in the 
denominator of Eq. (A9). As long as the real parts of the poles do not coalesce, i.e no 
two momenta li — lj and 1^ — lj are collinear, the angular integration still yields a factor of 
0(1) and the estimate (A10) remains correct. This is the regime that we are considering in 
the present appendix. The rest of the paper is devoted to the opposite limit, where all the 
external momenta are nearly collinear and the angular integral in Eq. (A9) is dominated by 
the vicinity of two antipodal points on the Fermi surface to which the external momenta are 
tangent. This observation motivates the introduction of the two patch theory in Section II 
and all the subsequent development of the present work. 

Returning to the non-collinear regime, upon combining Eq. (A10) with the Hertz scaling 
(A7), we conclude that the n-th term in the series (Al) scales as s n ^ 2 ~ l . Therefore, all terms 
with n > 2 represent non-local irrelevant perturbations, which confirms that the fluctuations 
with non-collinear momenta decouple. 

We would like to point out that the argument above still holds if one dresses the fermion 
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propagator by the one-loop self-energy, £(u;, k) ~ — isgncu|a;| 2//3 . This modifies the frequency 
dependence in the denominator of Eq. (A9) via, — i(l iT — lj T ) — > E(fc T + Z ir ) — T,(k T + Z jV ). 
However, since S(o;) <C wf|<t| f° r typical a; ~ |g| 3 , the estimate (A10) is still correct. 



Appendix B: Computations of Feynman diagrams 

Here we provide some details of the computations of the diagrams in Section V. 



1. Boson self-energy 

We begin by evaluating the two-loop polarization correction in Fig. 4, 

6 2 U(q) = Nj2 J <1 J0 ^f ) l 3 D(l)Gs(p)G s (p + q)G s (p - l)G s (p + q - I) (Bl) 

The contributions to the integral from the two patches are equal. Thus, integrating over p x , 
l x we obtain, 

dp T dp y dl T dl y 9(p T ) - 8(p T + q T ) 



5 2 U(q) — 2N [ d Prdpydl T dl 
W J (2vr) 2 (2tt) 2 



f({p} -{p + q}) + 2q yPy + q x + q 2 y 

x 6{l T -p T )-6(l T -p T -q T ) 

l jf({l - p - q} - {I - p}) + 2q y (p y - l y ) + q x + q 2 y 

where here and below we use the notation {p} = sgn(p r )|p r | 2 / 3 . We observe that the poles 
of the p y integral are always in the same half-plane. Thus, S 2 U(q) = 0. This is consistent 
with Ref. 50, which found that the two loop corrections to Eq. 3.3 are suppressed by factors 
of |w| 2 / 3 or M/lgJ ~ |cj| 2 / 3 . 

Now, let us proceed to compute the Aslamazov-Larkin diagrams, Fig. 6. We begin by 
evaluating the three point-function f s (q, I, —(l+q)) in Eq. (5.4). Note that f-(q, I, —(l+q)) = 
f + (P x q, P x l, —P x {l+q)) where P x (k , k x , k y ) = (k , —k x , k y ). The calculation of / is simplified 
when q T = 0. Then, performing the integral over p x and, subsequently, p y , in Eq. (5.4), 

f fq I -H + q )) qr =° f d Pr d Py *( (Pt + - 0(p T )) 



( 27r ) 2 ^({p + 0-M)-^-2/ y p y -/ 2 

i 

^{{P + 1} - {p}) -q x -l x - 2(q y + l y )p y + q*- l 2 y 
J_ f dp T \6(p T + l T ) - 6(p T )\(6(l y ) - 6(g y + l y )) 
2q y J 2tt zp. {{p + /} _ {p} ) + k _ + l y{qy + L 



(B3) 
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Thus, 

X + X-N 2 f dl T d 2 T dp T dp' T 



tSt-t/ „ -a A + A_iV 2 f dl T dH dprdp'.j, . 

^ = °< « = y + g) x 

|0(p T + Z r ) - <9(p T )||<9(p; + l T ) - 6{p' T )\Wy) ~ 0(ly + ft) I 



^({P + /} - {?}) + - ^Z y + ly(q y + ly) 



at({p' + - M) + f y k - Uft + k) 

+ (<?-> -g) (B4) 



^£({p/ + 1} _ {p/}) _ / x + + ^ + / 



Finally, integrating over Z a 



r2 



c--?^/ -j, A + A_A^ f dLdL dp T dp' ^ . ,. ^ . , 

* 3 n(ft = o, = -±^_ y -^££dw<! + fl ) x 



isgn(Z T )|%r + Z r ) - 0(Pr)|| W + Zr) " ^(Pv)ll^(O " ^ + ft) I >< 
1 



^({p + 0-W + {p' + /}-M) 

i 



-tc* 



+ (<? -> -?) 



at ({P + " M + {'P' + V- iP'}) + 2Z,(g, + l v ) 

(B5) 

The integral is invariant under q — > —q. Moreover, the integrals in the regions Z > and 
Zo < are related by complex conjugation. Thus, 

3 n ^ A+A_iV /-°°dZ T / ,,T dp T [ 1t dp' T f M dl y 1 1 



ft 2 7o 2vr J 2n J 2n J 2tt r Jr , c , ^ i X 
y y jo </o ./o \ + ? C6 |? s H ! , 1 ? 



C/ ((Z - P)? /3 + p' /3 + (Z - p')? /3 + /3 ) 

cjiij - pfj 3 + P 2/3 + (z - P o 2/3 + P ; 2/3 ) 

Cf({l - P)t 3 + Pr /3 + (Z - P')r /3 + Pr /3 ) 2 + *N*l*{\q y \ - l y f 



(B6) 



Notice that the integral over l y is bounded by the external momentum q y . This leads to a 
violation of the naive power counting, which would predict that each diagram in Fig. 6 has 
a superficial degree fo divergence A 2 ~ A 2 ^ 3 . Instead, we find that for Z r — > oo, the two 
diagrams behave as, 

/A \ 1/3 

5 3a U(0,q) = -5 3b n(0,g) ~ -X + X_N\q y \ -f (B7) 
(In reality, the divergence is cut once we exit the two patch regime where the momentum 
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l x <C ly. This occurs when l x ~ 3 becomes of order l y . However, for the Aslamazov- 
Larkin diagrams the internal momentum l y is controlled by external momentum q y . Hence, 
A r ~ q 3/2 and 5 3a U = -5 3b U ~ q 3/2 , as found in Ref. 26). 

However, as expected for problems involving a boson field coupled to the charge sector 
of the Fermi- surface, the divergence cancels when we add the two diagrams. In fact, for 
N ^> 1, the divergence is cut-off at jjlV 3 ~ q 2 , i.e. 



l T ~ N 3 ' 2 q 3 y /e 2 (B8) 

so that 2 

5 3 n(0,g) ~-A + A_iV 3 / 2 % (B9) 

e 2 

Note that the result is parameterically larger in the large-TV limit than the bare boson 
polarization, Eq. (2.7) (although it has the same scaling as the bare term). Also observe 
that the sign of the contribution (B9) is positive for the spin-liquid and negative for the 
Ising-nematic transition. 

One may ask whether the enhancement in (B9) is an artifact of taking q T = 0. However, 
since the integral in Eq. (B9) is saturated in the region (B8), we expect the result (B9) to 
be valid for, q T <C N 3 ^ 2 q 3 /e 2 , which is certainly satisfied by the typical bosonic momenta 
<?r ~ q y /e 2 . 

We can compute the proportionality factor in Eq. (B9) in the large- N limit. Changing 

/ \3/2 

variables to l T = Ijr) \q y \ 3 l T , p T = l T x, p' T = l T x', l y = \q y \y, 

Q 2 

5 3 n(0,g) = CA+A„-% (BIO) 

e 2 



c = - 25/W2 r * / W * /' *■ eV(1 - vY 



7T 



JO JO 



/ \ 3/2 _ / N 3/2 



X ^CR " ( B11 ) 

A(AH A J 3 + 4y 2 (l - y) 2 ) 

with, 

A = x 2 / 3 + (1 - x) 2 ' 3 + x' 2 ' 3 + (1 - x') 2 ' 3 (B12) 
For N ^> 1, the integral over l T is saturated in the region l T ~ 1, so, 



25/233/4^3/2 yoo ^ yl ^ yl ^ yl ^ y 3(! _ y ) 

Ti 



3 



C « / dx dx' dy - (B13) 

io ?J 3 Jo Jo Jo A(A 2 lt + % 2 (1 — |/) 2 ) 
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After a change of variables, z = ATJ 3 j{2y{\ — y)), 

(B14) 

The integral over x, x' can be evaluated numerically, 

, 1 



o jo 



dx I dx'j^ = 0.269653 (B15) 



so that 

C « -0.09601iV 3/2 , N ^roo (B16) 
We may also compute the constant C in Eq. (BIO) for the physical value N = 2, 

C « -0.04455 (B17) 



2. Fermion self-energy 



We next compute the three loop corrections to the fermion self-energy in diagrams Fig. 
10 b), c): 

5 36 E( Pr = 0,p) = iVA 3 A 3 J d j0± d j0 d jt0 G + ( P - h)G + (p - l 2 )G_{k)G_{k + h)G_(k + l 2 ) 
x DiljD^Dih-k) (B18) 

5 3c S( Pr = o,p) = iVA 3 + A 3 | + 'i)G+(p + / 2 )G„(A;)G4A: + h)G.(k + k) 
xD(h)D(l 2 )D(h-l2) (B19) 

Integrating over l\ x and Z 2 x we obtain, 

PYXpr = 0,0) = -iVA+A_ / ^ DWDtoDfr - l 2 ) A, 

y (2tt) 3 (2tt) 2 (27T) 2 _^^ /3 + <Jfc 
fl(ll T + fcj - gHlr) g(/ 2r + fc r ) ~ 0H 2r ) 

-f ((ii + kfj 3 + #*) + 5,7 + 2(k + p) y l ly - 6+ -^((/ 2 + kfj 3 + lT) + ^ + 2(* + - 5+ 

(B20) 
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8^(p T = 0,p) = -NX + X^ [ d ^ d -^%^D(l 1 )D(l 2 )D(l 1 - k) * - 

J (2tt) 3 (2tx) 2 (2tx) 2 _!££fc2/ 3 + 5 _ 

-f ((«! + A;)? /3 + ) + £ + 2(* +py iy + 2i?„ + 5+ 

9{l 2 r + fcr) ~ 0{-kr) 

-M 1 * + k)l /3 + it) + K + 2(* + p),/ 2 , + 2/ 2 2 y + 5+ 

(B21) 

where = ±p x + p 2 . Note the cancellation of the fermi-surface curvature terms l\ y ^ 2y in 
the "planar graph" 5 3b T,. 



We can reduce the integration range to k T > 0, as the region k T < is related by 
complex conjugation. There are then four different kinematic regimes: i) l\ T > 0, l 2r > 0, 
ii) l\ T < —k T , l 2r > 0, hi) l\ T > 0, l 2r < —k T , iv) l\ T < —k T , l 2r < —k T . The integral over k x 
in the regime i) vanishes as all the poles are in the same half-plane. The regimes ii) and iii) 
are related by l\ -h- l 2 . Thus, 



5 3b J:(p T = 0,p} = -iVA+A_ 



2 



oo 



dk T f d 2 k f°° dh T f°° dl 2r f dl ly dl 2y 



D(h)D(l 2 )D(l lT + l 2T J ly -l 2y ) — 

— jvfc- + S k 

1 



2tt J (2tt) 2 J kr 2tt J 2tt 7 (2tt) 
1 



^((/! - fc)^ /3 + # 3 ) + + 2(* + p)^i, - 5+ -^((/ 2 + ft)? 73 + # s ) + 5,7 + 2(* + p) y ky - 5+ 
dk T f d 2 k f°° dl\ T f°° dl 2r f dl\ y dl 2y 



l l 

!££ ((/l _ k) W + /J/3) + 5- + 2 (k + p) y l ly - 5+ l jf((l 2 - k) 2 T /3 + 1%') + 5- k + 2(fc + p),/ 2 , - fi+. 
+h.c. (B22) 
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S 3c E(p T = 0,p) = -N\ + \_ 



D(k)D(l 2 )D(l lT + hr, hy — hy)- 



dikj- f d k 



°° dl lT r °° 



o 2tt J (2vr) 2 J kr 2tt J 2tt J (2vr 



dX^r f dirndl 



1 



%({h - kfj 3 + /f ) + 5 k + 2(fc + p) y l ly + 2Z 2 y + 5+ 

1 

-W2 + k)l /3 + W) + £ + 2(fc + P )yhy + 2/ 2 2 y + 5+ 



dk T f cPk f°° dl\ T f°° dl 2r f dl\ydl 2 y 



2tt J (2vr) 2 J kr 2vr J kr 2vr 7 (2tt) 
1 



D{h)D{l 2 )D{h - l 2 ) 



1 



$((/ x - kfj 3 + /J/ 3 ) + 5~ + 2(k + p) y / lw + 2l\ y + *+ 

1 

l jf((h ~ W + llr) + *fc + 2(fc + V)yhy + 2^ + 5+ 



+ /i.e. 



N K T -t O k 



(B23) 



Integrating over k x and shifting k y ^ k y — p, 



S 3b E(p T = 0,p) = -A^A+A. 



2 i 



dfc T /" dfcy /"°° (i/i T /"°° g^t f dl\ydl 



2tx J 27r J kr 2tx J 2-k 



hr, hy — hy)' 



'2y 



(2^ 



(/St /3 + (h ~ k)l /3 + /if) - 2k y hy + 5+ 



+ 1 



-T$-((h ~ kfj 3 + l% 3 + (h + A;)? /3 + ill 3 ) + 2A^2 - Zx) w 
dhq- f dh, lt f dl^T f dlo-r- f dl-\*,dl 



2tt / 2ir L 2-k J k 2vr 



^D(h)D(l 2 )D(h-l 2 ) 

1 



^(k^ 3 + (I, ~ kfj' + ll' 3 ) + 2kyl ly - 6+ ^(0 3 + (1, - kfj' + IT ) + 2^/2, " S p 



h.c. 



(B24) 
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5 3c £(p T = 0,p) 



-iVA+A- 



2% 



2ix J 2ix J kr 2ir 



10 dl\ T f°° dl 2r f dh y dl 2 y 



2tt 



(27T) 



D(h)D(l 2 )D(h T + hr, hy — hy)- 



1 



icf,, 2/3 



JV -(r + (li-fc)f + ff)-2V 
1 



+ i 



mh - k) i j A + + (/ 2 + k)r + + 2^(/ 2 - h) y + 2(/ 2 \ - 1%) 

dl^T (* dirndl*) 
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oo 



dkt T f dky j" dl\q- 



2tt ./ 2tt ./fc T 2tt y fcT 2vr 



(27T) 



D(h)D(l 2 )D(h-l 2 ) 



^(0 3 + (h - kfj 3 + l% 3 ) + + 2P ly + 5+ ^ /3 + (l 2 - kfj 3 + if 3 

+h.c. 



2kyl 2y + 2ll y +S+ J 



(B25) 



The integration regions l\ y > and h y < give the same contribution. So, integrating over 
k,, 



S 3b E(p T = 0,p) = iVA+A. 



D(h)D(l 2 )D(l lT + l 2T ,l ly - ky) 



dk 1 



d\\ T f dl 2r f dl-\i, f dl 



ly I u 'l 2 y 

lly 



+ 



-%{h y {{h - k)l' 3 + it 3 + 3 ) + l ly ((l 2 + kfj 3 + fj 3 - 3 )) + (l 2 - h) y S+ 
dk T r°° dl w f°° dk 
27 



~"ir f°° dl 2r f°° dl\y r°° dl 2y 



2tt 1 1. 2tx 



2ix J 2ix 
1 



D(h)D(l 2 )D(l lT -l 2T ,hy + hy) 



-t(hy((h ~ kfj 3 + if + 3 ) + l ly ((l 2 - kfj 3 + fj 3 + 3 )) + (h + l 2 )y8+\ 



+ h. 



(B26) 



5 3c Z(p T = 0,p) = iVA+A_ 



dk T f 00 dl lT f°° dl 2r f°°dl 



2ir J kr 2ix J 2tx J 2tt J hy 2ir 



'-ly f°° dl 2y 



D(h)D(l 2 )D(l lT + 

hri hy — hy) 



l jf(hy((h - k)T + ll /3 + kl' 3 ) + hy((h + k)T + IT ~ kl' 3 )) + 2l ly hy(h - h)y ~ (h ~ h)yS+ 



+ 



dkr f°° dl lT 



2tt J kr 2ir J kr 2tt 



dl 2 r f dl\y f dl 



2tt 



22/ 



2tt 



D(h)D(h)D(l lT - hr, hy + hy) 



1 



-^-(h y ((h - kyr + it + k' 6 ) + hy((h - k)r + i;i 6 + kr)) - 2i ly i 2 y(h + i 2 ) y - (h + i 2 ) y s+i 

+ h.c. (B27) 



Expanding the self-energy in 5+ and performing a change of variables h T = k T x±, l 2r 
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k T x 2 , h y = (c b e 2 k T ) 1 / 3 y 1 , l 2y = (c b e 2 k T ) 1 / 3 y 2 , 

r°° dk 

5 3b Z + (p T = 0,p) = A + A_(Jx + J 2 )5+ / — ^ (B28) 

JO fc T 

5 3c S + (p T = 0,p) = 5 3c Z + (p T = 0,p = 0) + A+A_(J 3 + J4)S; [ 

Jo 



T 1 (B29) 

IV n- 



where 



Ji = f dx x [ dx 2 [ dyi [ dy 2 - 



7T 2 Ji Jo Jo J yi {xi + y\){x 2 + yl){xi + x 2 + {y 2 - yi) 3 ) 

x ; 27s ; 27s ( B3 °) 

(y 2 ((x 1 - l) 2 / 3 + xT + 1) + Vl ((x 2 + l) 2 / 3 + x 2/3 - l)) 2 

J 2 = ^- [ dxi [ dx 2 [ dyi [ dy 2 - 



7T 2 Ji Ji Jo Jo Oi +Vi)( x 2 + y 2 )(\xi -x 2 \ + ( yi +y 2 ) 3 ) 

x ; 27s ; 27s ( B31 ) 

(y 2 (( Xl - l) 2 / 3 + xT + 1) + yi((x 2 - I) 2 / 3 + xT + I)) 2 
6 r*> [~ f° [°° ym(y2-yi) 2 

^3 = 37775- / dxi I dx 2 I dyi \ dy 2 — — —3-- — 3 



tt 2 N 2 J 1 J J J yi ( Xl + yf)(x 2 + y|)(xi + x 2 + (y 2 - 7/1) 3 ) 

x HvKyi - yi? - Mv^i - 1) 2/3 + 4 /3 + 1) + yi((x 2 + 1) 2 / 3 + xf - 1)) 2 

* (3y 2 y 2 (y 2 - Vl ) 2 + ^(y 2 {(x x - If I 3 + xf + 1) + yi ((x a + l) 2 / 3 + x 2/3 - l)) 2 ) 2 

(B32) 

J a = — 7-777 / dxi / cfe 2 / cfo/i f dy 2 



7r 2 iV 2 7i ' 7i Jo Jo (x 1 +yf)(x 2 + yl)(\x 1 -x 2 \ + (yi+y 2 ) 3 ) 

x + y 2 f - ^M(xi ~ 1) 2/3 + + 1) + yi((x 2 ~ 1) 2/3 + + l)) 2 

X {Zylvliyx + y 2 ) 2 + Mv^i - i) 2/3 + 4 /3 + 1) + lfr((* 2 - 1) 2 / 3 + x 2 2 /3 + l)) 2 ) 2 

(B33) 

Cutting off the UV divergence in (B28), (B29) at k T = A T ~ A 3 /e 2 , we obtain to 
logarithmic accuracy, 

S 3b J: + (p T = 0,p) = A + A„(J! + J 2 )5+log-^ (B34) 

A 3 

5 3c T, + {p T = 0,p) = 5 3c T, + {p T = 0,p = 0) + A + A_(J 3 + JJS} log (B35) 

which is equivalent to Eqs. (5.9), (5.10) with J b = 3(Ji + J 2 ), J c = 3(J 3 + J4). Note that Ji 
and J 2 are constants independent of N, 

Ji « 0.01276 (B36) 
J 2 » 0.02264 (B37) 
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On the other hand, the constants J3 and J4 are N dependent. In the large- N limit we can 
evaluate these constants analytically to leading logarithmic accuracy by setting N = 00 in 
the integrand. 

2 /*oo roo poo /*oo ^ 

1/3 ~ n 2 N 2 J 1 dXl J dX2 J dyi J yi ^ymfa + yf)(x 2 + + x 2 + (y 2 - Vi) 3 ) 

(B38) 

The above integral diverges logarithmically when yi,y 2 ,x 2 — > 0. Hence, 



— t- / dx 2 dy x \ dy 



tt 2 N 2 J 1 x\ J """'Jo J yi ^ y2 yiy 2 (x 2 + y%) 



2 f 1 dy 2 , , _ 3 , f V2 dyx 



dx x f 1 . f 1 . f 1 , 1 2 f 1 dy 2 

■^E{V 2 ') I 
Jo 

(B39) 

Inspecting the original integral (B32), we observe that the logarithmic divergence in (B39) 
is cut-off when yi(y 2 — yi) ~ j*. Hence, 



n 2 N 2 J N -\ y 2 J(Ny 2 )~i Vi 4vr 2 iV 2 



Similarly, 



/*oo /*oo /*oo /*oo 

rfXx / cfe 2 / d Vl I dy 2 



n 2 N 2 J 1 7i Jo Jo Vm{xi + yr)(^2 + vl){\xi - x 2 \ + {y x + y 2 f) 

log((yi + y 2 )- 3 ) (B41) 



4 pdx! f 1 ^ pdyi,. 



7r 2 iV 2 J 1 x\ Jo y 2 Jo V\ 

Inspecting Eq. (B33), we see that the logarithmic divergence in (B41) is cut-off when 
2/12/2 ~ jf- Writing, y x = y 2 z, 

12 f 1 dy 2 f 1 



n 2 N 2 



I dm 1 dz + + ^ i^ log3iV (B42) 



We note that expressions (B40), (B42) do not include subleading polynomial corrections in 
logiV. We can also calculate the constants J3, J4 numerically for N = 2, 

J 3 ~ -0.004491 (B43) 
J 4 » -0.008158 (B44) 



Finally, we compute the insertion of the 2 operator into the fermion two-point function, 
which determines the renormalization of the chemical potential 5 away from criticality. The 
UV contribution at three loop order comes from the diagrams in Figs. 10 b) c) and can be 
obtained by expanding the bosonic propagators in Eqs. (B26), (B27) to linear order in r. 
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This yields, 



5 s b 9Z uv _ N 



Or 



dk T 



a dh T f°° dh T f°° dl\y 



dl 



2y 



'0 Jk T * n JO * n JO * n Jhy 

{D{l l ) + D{l 2 ) + D{l lT + 

hri hy ~~ 

l 2y ))D(h)D(l 2 )D(l lT + l 2r ,l ly - hy) 

1 



ic 



N Hky((h - k)T + iT + tr) + i ly ((i 2 + k)r + it - kr)) 

00 dk T f°° dhr f°° dh 



,2/3> 



>2/3 , ,2/3 



,2/3n 



+ 



""lr /"°° <^2r /"°° ^lj/ /"°° dl 2y 



27T ,/ feT 27r J kr 2tt J 2tt Jo 2tt 
(D(k) + D(l 2 )+D(l lT - 

hri hy 

+ l 2y ))D(l 1 )D(l 2 )D(l lT -l 2r ,l ly + hy) 

1 



-^f(M(/i - ^)? /3 + i% z + ^ /3 ) + i ly ((h - k)r + ^ + kY A )) 



,2/3 



2/3 ,2/3 



,2/3n 



+ /l.C. 



(B45) 



°° C?Zl r (iZ 2l - /" f dl 2 y 



(D(k) + D(l 2 ) + D(l lT + 

hri hy ~ h y ))D{h)D(l 2 )D(hr + h T ,hy - hy) 
1 



f (/*((«! - A:)? /3 + # s + ^ /3 ) + l ly ((h + ^)r /3 + if - kt)) + 2/ ly / 2 ,(/ 2 - h) 
dl 2 T f°° dhy f°° dl 2 y 



dk T [°° dl lT r °° 



2vr j kr 2vr J kr 2tt J 2tt ,/„ 2vr 



(B46) 



(D(h) + L>(Z 2 ) + £>(Z lr - l 2 r, hy + hy))D(h)D(l 2 )D(hr - hr, hy + hy) 

1 



■^(h y ((h - kfj 3 + if + 3 ) + h y ((h - kfj 3 + I 2 ' 3 + 3 )) - 2Z 1 ,Z 2 ,(/ 1 + / 2 ) 



+ h.c. 



(B47) 



We observe that the contribution from the diagram in Fig. 10 b) vanishes, while the diagram 
in Fig. 10 c) gives upon switching to dimensionless variables, 



o<9S uv T o 1 a 
dF = g s 



(B48) 



13 



with 

. 36 y°V r r°° y~ y^f - yi) 2 

/ yi , 1/2 _^ z/2 - yi \ 

\xx + yf x 2 + yl x 1 + x 2 + {y 2 - yi) 3 ) 

1 

32/?2/I(2/ 2 - yi) 2 + ^(y 2 ((*i - 1) 2/3 + x? /3 + 1) + yi((x 2 + l) 2 / 3 + xf - l)) 2 

18 r, ^Kyi+1/2) 2 

( Vl J_ ^ 2 , 2/1 + 2/2 \ 

Vxi+y 3 x 2 + yl \xi - x 2 \ + (y 1 + y 2 ) 3 J 

1 

3i/! 2 y 2 2 (yi + y 2 ) 2 + ifcWfa - i) 2/3 + xT + l) + yi((x 2 - i) 2/3 + * 2 /3 + 1)) 2 

(B49) 

Evaluating the above integral, we obtain Eq. (5.18). 
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